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I. INTRODUCTION 

For more than 30 years the microscopic descrip- 
tion of nuclear ground state properties is relying on 
the self-consistent mean field approach, or Hartree- 
Fock (HF) approach built from effective in-medium 
nucleon-nucleon interactions Q. The mostpopular 
interaction is the Skyrme-type interaction^ whose 
analytic form leads to considerable simplifications 
of the HF calculations in finite nuclei. The gen- 
eral point of view is that the Skyrme interaction 
is a phenomenological one whose parameters are di- 
rectly adjusted on a few selected observables taken 
from infinite matter and some doubly-magic nuclei. 
In its most sophisticated versions, the Skyrme-type 
force can predict binding energies of all measured 
nuclei with an overall error of less than 0.7 MeV[3- 
Other parametrizations try to improve the descrip- 
tion of systems with a large neutron excess by 
incorporating constraints from variational calcula- 
tions performed for neutron-rich and pure neutron 
matter IT 



On the other hand, important progress have been 
made in the recent years concerning Brueckner- 
Hartree-Fock (BHF) calculations of infinite matter^, 
8]. For many years, a major drawback of the BHF 
approach was that the empirical saturation point of 
symmetric nuclear matter could not be reproduced 
if one starts with only a realistic two-body bare in- 
teraction (the so-called Coester line problem). This 
is the reason why in the early attempts to derive 
from a Brueckner G-matrix an effective interaction 
suitable for HF calculations of finite nuclei it was 
necessary to renormalize phenomenologically the G- 
matrix[9lllC|. Two recent achievements have led the 
BHF approach to a quite satisfactory status: one is 
that the convergence of the hole-line expansion has 
been proved to occur already at the level of three- 
body correlations 7J; the other one is that the inclu- 



sion of three-body forces improves to a large extent 
the prediction of the saturation point 8] . 



Thus, it is timely to re-examine how one can re- 
late the BHF description of homogeneous matter to 
a HF description of finite nuclei. Our strategy is to 
look for a Skyrme-like parametrization adjusted so 
as to reproduce the BHF results calculated in sym- 
metric nuclear matter and also spin- and isospin po- 
larized matter, then to examine the HF predictions 
of this parametrization in nuclei. In our approach 
it is possible to determine all the parameters of the 
force except the two-body spin-orbit parameter Wo 
because this component of the force gives no contri- 
bution in homogeneous matter. This will remain an 
adjustable parameter when calculating finite nuclei, 
and its value is fixed on the lpl/2 — lp3/2 splitting 
in 16 as it is usually done. Of course, the two-body 
Coulomb interaction has also to be added. 



In a recent work, Baldo et al. 11] have studied 
the relation between the Brueckner results in infi- 
nite matter and some of the Skyrme parameters, but 
the number of constraints used was insufficient to de- 
termine the velocity-dependent part of the Skyrme 
force, which is a very important part since it governs 
the effective mass behaviour. 



In this work, we start from the recent BHF calcula- 
tions of Refs. [a, uM extended to spin- and isospin po- 
larized homogeneous matter. The paper is organized 
as follows: in Sec. II a brief review of relevant BHF 
predictions is given. In Sec. Ill we explain how the 
Skyrme parameters are determined. Sec. IV presents 
and discusses the HF results obtained for selected 
closed-shell nuclei. Concluding remarks are given in 
Sec.V. 
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II. REVIEW OF THE RELEVANT BHF 
PREDICTIONS 

For densities around the saturation value and 
smaller, the BHF approximation, in which the en- 
ergy shift is truncated at two hole-line level, is a quite 
good approximation. In that region in fact we are al- 
lowed to neglect the three-body correlations, which 
anyhow give a small contribution when adopting the 
continuous choice for the auxiliary potentialQ. On 
the other hand, the BHF approximation embodies 
both the inert core pure BHF mean field and the core 
polarization term, the latter arising as rearrangement 
potential^]- Thus, the main effects of the correla- 
tions on the effective mass will be taken into account. 
The fit of the microscopic effective mass is a qualify- 
ing aspect of the present parametrization, since this 
quantity is one of the least well constrained proper- 
ties of the phenomenological Skyrme forces |13j. 

Let us briefly review the salient features of the 
BHF equation of state (EoS): 

i) saturation point 

The BHF approach with only two-body force size- 
ably overestimates the saturation density^]- This 
is commonly attributed to the missing effect of the 
three body force in the region above the empiri- 
cal saturation density. After introducing a three 
body force in the calculations the saturation den- 
sity is improved without appreciable change of the 
corresponding energy. In the most recent Brueckner 
calculations j&j it was estimated , po ~ O.I8/to~ 3 , 
still beyond the range of the empirical values of 
the central density of nuclei, and the energy -j ss 
— lA.8MeV. The main effect of such a failure is ex- 
pected to appear in the calculation of the neutron 
and proton density profiles. 
ii) effective mass 

The nucleon effective mass inside the medium is an 
outcome of the BHF approach and as such, it must be 
reproduced by the equivalent Skyrme parametriza- 
tion. It is related to the momentum dependence 
of the on-shell self-energy E(efe,fc), where €k is the 
quasiparticle energy. It plays an important role not 
only in the theoretical description of the transport 
phenomena, including heavy ion collision (HIC) sim- 
ulations, but also for the level density of nuclei. Re- 
cently the isospin splitting of the effective mass has 
been the subject of some debate since different ap- 
proaches predict contradictory results [l5l[l^ . [r^ . [lq . 
In the framework of the Brueckner approach it is 
possible to trace the isospin splitting of the effec- 
tive mass back to properties of the bare nuclear in- 
teraction 0, ^| . It turns out that the neutron ef- 
fective mass is linearly increasing with the neutron 
excess while the proton effective mass is symmetri- 
cally decreasing. Embodying the neutron and proton 
effective masses in the fit one may expect important 
isospin effects in observables which are sensitive to 



the effective mass itself. 



Hi) symmetry energy 



The symmetry energy a s (p) has stimulated a lot 
of interest for its relevance in HIC physics, nuclear 
astrophysics and exotic nuclei. In fact it is related 
to the isospin splitting of effective mass, neutron and 
proton mean fields, etc... In particular, the neutron 
skin in neutron-rich nuclei seems to be very sensitive 
to the details of the density dependence of a s (p) [23] • 
In the microscopic approaches, the saturation point 
becomes less and less stable with increasing neutron 
excess and at some critical point, before reaching 
the conditions of pure neutron matter, it disappears. 
This transition formally amounts to the transition 
from a minimum to an inflexion point in the func- 
tion Correspondingly the symmetry energy 
would also exhibit an inflexion point as a function 
of density. This entails that any parametrization of 
a s (p) such as p a , which is often adopted in calcula- 
tions, is not suitable to reproduce this behavior. This 
seems to be the case with the Skyrme forces. In the 
BHF approximation the symmetry energy at the sat- 
uration point turns out to be a s (p ) « 34.MeV Q. 
At low density it is independent of the force and 
therefore it can be considered well established from a 
quantum-mechanical many body theory. Above the 
saturation point the three body force has a strong 
influence on the symmetry energy. The BHF predic- 
tion is in a rather good agreement with the relativis- 
tic Dirac-Bruecknerll (J but both diverge from vari- 
ational calculations 5, 6] and also from some Skyrme 
forces (for a discussion, see Ref.[r3jp. The structure 
of neutron stars has been addressed as a possible con- 
straint for the EoS of asymmetric nuclear matter and, 
in particular, for the symmetry energy, but so far the 
calculations do not give a definite answer, 
iv) Landau parameters 

In order to have a full determination of the Skyrme- 
force parameters and not only some combinations 
of them, it is not sufficient to fit bulk properties 
such as binding energies, effective masses and sym- 
metry energy, as we shall see in Sec. III. We will 
use the additional constraint of reproducing the Go 
Landau parameter extracted by extending the BHF 
calculations to spin and isospin asymmetric nuclear 
matter j2JJ. The Brueckner predictions for the Lan- 
dau parameters plj have proved to reproduce the ex- 
isting experimental data; in particular the parame- 
ter G' is consistent with the centroid energy of the 
Gamow- Teller giant resonance j2^. So far only the 
values of the Landau parameters at the saturation 
point can be tested. At lower density experimental 
information can come from the study of giant reso- 
nances in exotic nuclei. 
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III. DETERMINATION OF SKYRME 
FORCE PARAMETERS 



obtained could describe reasonably well finite nuclei. 



The standard form of the Skyrme effective inter- 
action is: 

7(ri,r 2 ) = *o(l + x P a )8{v) 



1 



h{\ + x 1 P a )[P' 2 S(r) + S{r)P 2 } 



+t 2 {\ + x 2 P a )V ■ 8{r)V 
+iW a • [P' x S(r)P] , 



(1) 



where r = ri — r 2 is the relative coordinate of the 
two particles and R = (ri + r2)/2 the center of mass 
coordinate. P = (Vi — V 2 )/2i is the relative momen- 
tum acting on the right, and P' its conjugate acting 
on the left. P a = (1 + o\ ■ ai)/2 is the spin-exchange 
operator. 

The force parameters are the ti and Xj, the power 
a of the density dependence, and the spin-orbit 
strength Wq. As already mentioned, the latter pa- 
rameter will be adjusted phenomenologically in some 
specific nucleus. As for the a parameter, it is diffi- 
cult to extract it in a unique way just from nuclear 
matter bulk properties. In the literature there are 
several classes of Skyrme forces characterized by the 
value of a, the most common values being 1/6, 1/3, 
1. Since we do not have enough constraints from 
BHF calculations of nuclear matter to determine all 
the free parameters, we choose to adopt in this work 
cr=l/6 and to concentrate on finding the remaining 
parameters. 

The basic inputs of this work are the results of 
the BHF self-energy which includes the core polariza- 
tion term (called in the literature extended BHF ap- 
proximation (EBHF)[j3) and the EOS with two and 
three-body forces^ . The procedure to determine the 
force parameters proceeds through three main steps. 
The first step concerns the fit of the nucleon effective 
mass in symmetric and non-symmetric nuclear mat- 
ter. This enables one to find the values of two im- 
portant combinations - O s and 0„ - of tx, fa, X\, x^. 
In the second step we look at the energy per par- 
ticle in symmetric and non-symmetric nuclear mat- 
ter as a function of total density and neutron-proton 
asymmetry. The fit of these quantities determines 
several families of (to, ts, xq ,Xs) parameters. In the 
last step, we make use of the constraints imposed 
by the value of the Go Landau parameter and by 
a combination of parameters governing the surface 
properties of finite nuclei. In this way, the remain- 
ing parameters t\,t2,Xi, x% are uniquely determined 
for each parameter family since we already know the 
values of S and <d v . Because the fit of BHF bulk 
properties has been supplemented with the surface 
condition one may hope that the parameter sets thus 



A. Effective masses 

In symmetric nuclear matter, the isoscalar effec- 
tive mass of a nucleon has the following expression 
(here and in the rest of this paper we follow the same 
notations as in Ref.0): 



m; . m ^ 

where Q s = [3i a + (5 + 4x 2 )t 2 ] 
an isovector effective mass as 



(2) 



One can also define 



ml , m „ , 

— = 1 + TZipBv y 
m 4n z 



(3) 



where Q v — ti(xi + 2) + t 2 (x 2 + 2). In asymmetric 
nuclear matter with an asymmetry parameter (3 = 
(N — Z)/A, the nucleon effective mass is 



= i 1 + ^5P @ s ~ Z^^Qv - Q s )Pp) 



8h 2 



8h 2 



(4) 



with q = 1 for neutrons and q = — 1 for protons. 

To obtain Q s we fit the values of the nucleon effec- 
tive mass calculated in symmetric nuclear matter in 
EBHF approximation [12 with the three-body force 
effects included. The fit is illustrated in Fig.l, and 
the resulting value from the fit is <d s = 400.8 MeV 
fm 5 . Next, we can determine Q v by fitting m*/m cal- 
culated in asymmetric nuclear matter. The fits are 
shown in Fig. 2, and the corresponding value of 0„ is 
356.4 MeV fm 5 . Actually, the fits of m*/m remain 
good beyond (3=0A . 
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FIG. 1: Nucleon effective mass m*/m in symmetric nu- 
clear matter. 

At this point we can already make an interesting 
observation. Some Skyrme parametrizations of the 
literature predict that, for a fixed density, the pro- 
ton and neutron effective masses are respectively an 
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increasing and decreasing function of the asymme- 
try parameter /J^jjJ. This behaviour is opposite to 
that predicted by EBHF calculations |15| . In Fig. 3 
we show a comparison of effective masses calculated 
with the parameter set SLy4[4| and with our values 
of O s and 0„. Thus, our Skyrme parametrization 
will differ from some usual parametrizations as far 
as the isospin splitting of effective masses is con- 
cerned. Recently, it was pointed out0 that ef- 
fective masses obtained in Dirac-Brueckner-Hartree- 
Fock have a similar behaviour to that of EBHF, and 
opposite to that of the relativistic mean field (RMF) 
model. It must also be noted that, if one includes the 
Fock terms in a relativistic Hartree-Fock description, 
the trend of the neutron-proton mass splitting be- 
comes closer to that of the Dirac-Brueckner-Hartree- 
Fock and therefore, the behaviour of the RMF mass 
is due to the omission of exchange terms[25j. 
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FIG. 3: Neutron (upward triangles) and proton (down- 
ward triangles) effective masses m*/m calculated with 
SLy4 (in white) and with the present parametrization (in 
black) . 



B. Energy per particle 



In symmetric matter the energy per particle has 
the following expression: 
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FIG. 2: Nucleon effective masses m*/m in asymmetric 
nuclear matter, at two different asymmetries. Triangles 
are EBHF results, the lines are the fits. 
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,(5) 

More generally, in asymmetric matter characterized 
by an asymmetry parameter (3 the energy per particle 
is: 



3fa? , 3tt^ . a 2 „ 1 , . rtf)1 

3 P 3F f + -top 3 - (2x + I)/? 2 ] 

+ ^3P ff+1 [3-(2x 3 + l)/3 2 ] 



3 37T 2 2 5 

40 V 2 ; F 1 3 



2e v )Fs 



(6) 



where F m {/3) = |[(1 + P) m + (1 - /3) m ]. 

We first start with the case of symmetric matter. 
The fit of the BHF values of j;(ft = 0, p) determines 
several possible values for the couple (to,t 3 ). To limit 
the number of possible couples we impose that the 
compression modulus Koo must be between 200 and 
240 MeV, and the pressure P=0 at the saturation 
density pq. The optimal parameter sets we find have 
Koo around 210 MeV. As for p it is about 0.18/m~ 3 
in the BHF calculation with three-body force[j|, a 
value somewhat larger than that usually adopted for 
Skyrme parametrizations. The consequence will be a 
general underestimate of radii in finite nuclei as we 
shall see in Sec. IV . 

For each (to,ts) couple previously determined we 
add the constraints of fitting ^ 0,p) as well as 
the symmetry energy a s (p) of infinite matter calcu- 
lated in BHF. In terms of Skyrme parameters the 
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symmetry energy a s (p) is: 



1 8 2 {E/A) 

2 dft 2 



1/3=0 
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1 ft 2 ,371-2 2 . 

= 32^ ( — 8*0(2x0 + 1../^ 

1 . 37T 2 , 2 , _ _ . 5 

-^( — ) 3 (3e t ,-29 s ) P 3 



1 

48 



(7) 



In practice, we choose the EOS with ft = 0.4 as 
the fitting object. In the fitting procedure we set 
the symmetry energy at saturation density, a s (po) = 
34MeV as a constraint for the full EoS of asymmet- 
ric nuclear matter. Thus, we obtain several possible 
solutions for (to, Xo, £3, X3) with to, £3 already given 
by the symmetric matter fits and xo,x 3 determined 
by non-symmetric matter properties. Figs. 4-5 dis- 
play typical fits obtained in this way. In Fig. 4 are 
shown the energies per particle in symmetric and 
ft=0A non-symmetric nuclear matter, calculated in 
BHF and by the Skyrme procedure. The pure neu- 
tron matter case, which is not included in the fit- 
ting procedure, is also shown to demonstrate that the 
present Skyrme parametrizations describe reasonably 
well the BHF energies per particle in the whole range 
of ft from to 1 and for densities up to the highest 
BHF calculated values around 0.35 /to -3 . However, 
one can see that the Skyrme energy functional is not 
able to reproduce accurately the EoS of pure neutron 
matter which exhibits an s-shaped behaviour. This 
suggests that a change of the analytical structure of 
the energy functional would be necessary. In Fig. 5 
is shown the fit of BHF symmetry energy. The devi- 
ations observed in the neutron matter case in Fig. 4 
reflect in the fit of the symmetry energy since the 
symmetry energy is calculated from the expression 

a s{p) = JiipiP = 1) — §(Pj — 0)i an d then a bad 
fit for the neutron matter EOS induces a bad fit for 
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FIG. 4: Energy per particle as a function of density. 
Squares, dots and triangles are BHF results for symmet- 
ric matter, (3=0.4 asymmetric matter and neutron mat- 
ter, respectively. Solid and dashed lines are results of the 
fit, the dotted line is the extrapolation at /3=1. 
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FIG. 5: The symmetry energy from BHF calculations 
(squares) and the present fit (solid curve). 



To summarize, the bulk properties of nuclear mat- 
ter have enabled us to determine several sets of 
(to, xq, ts, X3) parameters and the values of the pa- 
rameter combinations O s and 0„ for a fixed value 
cr=l/6 of the density dependence. To complete the 
task we proceed to determine the remaining param- 
eters (ti, Xi,t2, X2) in the next subsection. 



Constraints from Landau parameters and 
finite nuclei 



The bulk nuclear matter properties give only the 
two parameter combinations O s and Q v . To proceed 
further, we can use some additional constraints from 



the Landau parameters corresponding to our refer- 
ence EBHF calculation. These Landau parameters 
have been investigated in Ref . [2j . The two parame- 
ters Fo and Fq are related to the compression mod- 
ulus and symmetry energy, respectively, which are 
already incorporated into the preceding constraints. 
Then, only Go and G' are left as constraints, and 
in principle they are enough to determine the four 
remaining parameters. However, we prefer to use 
as constraint only one of them and keep some free- 
dom to better optimize the surface properties of finite 
nuclei. Indeed, the surface effects cannot be deter- 
mined from infinite matter calculations and we need 
to adjust these surface effects by performing Skyrme- 
Hartree-Fock (SHF) studies of some selected nuclei. 
We will choose to fit Go whose value in EBHF|2lT| 
is 0.83 at po=0-18 /to -3 . The parameter Go is ex- 



G 



TABLE I: The Skyrme parameter set and the correspond- 
ing bulk properties of infinite nuclear matter. 

LNS 

t (MeVfm 3 ) -2484.97 

t^MeVfm 5 ) 266.735 

t 2 (MeVfm 5 ) -337.135 

t 3 (MeV/m 3+3<T ) 14588.2 

x 0.06277 

xi 0.65845 

x 2 -0.95382 

x 3 -0.03413 

a 0.16667 
W (MeVfm 5 ) 96.00 
p (/m~ 3 ) 0.1746 
E/A(MeV) -15.32 
K^iMeV) 210.85 
^-(isoscalar) 0.825 
^(isovector) 0.727 
a s (MeV) 33.4 



pressed in the Skyrme parametrization as: 

G = N Q {2A + 2{3ir 2 p/2)^B) (8) 

where Nq is the level density at the Fermi surface 
and 

B = -\tx(\-x 1 ) + ±t 2 {±+x % ). (10) 

In the SHF energy functional an important term 
governing surface properties of N=Z systems is the 
term0| a s X7 2 p where a s = (t 2 (5 + 4a;2) -9ii)/32. To 
determine the 4 parameters (t\, X\,t2, X2) satisfying 
the fixed values of O s and <d v , and also the spin-orbit 
parameter Wo we adopt the following procedure. We 
choose as reference nuclei the closed-shell and closed- 
subshcll nuclei 16 0, 40 > 48 Ca, 56 > 78 Ni, 90 Zr, 100 ' 132 Sn, 
208 Pb. We vary a s in a range of values similar to 
that of usual Skyrme forces, and for each value we 
obtain a set of (ti, x%, i 2 , #2) with which we can per- 
form a SHF calculation of the reference nuclei. The 
corresponding value of Wq is obtained by adjusting 
the pl/2 — pi/2 proton splitting in 16 0. In this way, 
we can determine the set which gives the best overall 
results for binding energies and radii in the reference 
nuclei. 

Table I summarizes the outcome of the fit. The 
full parameter set is called LNS. In the lower part of 
Table I are shown the main bulk properties of nuclear 
matter calculated with LNS. It can be noted that po 
is larger and the saturation energy is slightly less 
negative than the empirical saturation point (0.16 



/to -3 , -16. MeV). The consequence is that, in finite 
nuclei central densities will tend to be too large and 
radii become systematically underestimated. 



IV. HF CALCULATIONS OF MAGIC 
NUCLEI 

We now discuss the results of HF calculations of fi- 
nite nuclei made with the LNS parametrization. The 
parameter set of Table I is supplemented with the 
two-body Coulomb force. The Coulomb exchange 
contributions are treated in the Slater approxima- 
tion. The center of mass correction to the total en- 
ergy is approximated in the standard way, kee ping 
the one-body and dropping the two-body terms |26j. 
The HF equations are solved in the radial coordinate 
space, assuming spherical symmetry. 

We choose to test the LNS parametrization on the 
following set of closed-shell and closed-subshell nu- 
clei: 16 0, 40 Ca, 48 Ca, 56 Ni, 78 Ni, 90 Zr, 100 Sn, 132 Sn, 
and 208 Pb. The LNS force has not been fitted on 
finite nuclei and therefore, one cannot expect a good 
quantitative description at the same level as purely 
phenomenological Skyrme forces. It is nevertheless 
interesting to compare its predictions with those of 
a commonly used force like the SLy4 interaction. In 
Fig. 6 are shown the relative deviations of charge radii 
(upper panel) and energies per particle (lower panel) 
calculated with LNS and SLy4. As for the binding 
energies, one can see that LNS is doing reasonably 
well in the Ca and Ni regions but it is underbind- 
ing somewhat 16 O and overbinding the medium and 
heavy nuclei, the discrepancies remaining within a 
5% limit. The SLy4 force is doing of course much 
better since it is adjusted on doubly magic nuclei. 
The deviations of the LNS energies can be attributed 
to the incorrect saturation point of the BHF equation 
of state, and also to the lack of information concern- 
ing surface properties that one should fulfill. The 
charge radii of LNS exhibit a systematic behaviour 
of underestimating the data by 2 to 4%. This can 
be understood again as a consequence of the BHF 
saturation point being shifted towards a larger den- 
sity. Then, the central density in nuclei calculated 
with LNS becomes larger than what it should be and 
this reduces the spatial extension of nuclear densi- 
ties. This is illustrated in Fig. 7 where we show the 
neutron and proton distributions in 208 Pb calculated 
with LNS and with SLy4. 

Finally, we would like to comment on the spin- 
orbit component of the LNS parametrization since 
this is the only parameter that we could not relate 
to the EBHF calculation. We have fitted it to re- 
produce the experimental lpl/2 — lp3/2 spin-orbit 
splitting of neutron and proton levels in 16 O. The 
value W =96.MeVfm 5 that we find seems somewhat 
smaller than for other Skyrme forces where Wq usu- 
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cles, respectively. 
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FIG. 7: Neutron and proton densities in 208 Pb, calculated 
with LNS and SLy4 parametrizations. 
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FIG. 8: The spin-orbit potentials in le O. 



ally ranges from 105. to 130. MeVfm 5 . In Fig. 8 we 
show that the spin-orbit potentials in 16 O calculated 
with LNS and SLy4 are nevertheless very close and 
they must give the same spin-orbit splitting. 



V. CONCLUSION 

In this work we have derived a Skyrme-type 
parametrization of an effective interaction suitable 
for Hartree-Fock calculations of finite as well as in- 
finite systems. The starting point is the BHF cal- 
culations of infinite nuclear matter at different den- 
sities and neutron- proton asymmetries. These BHF 
studies include effects of three-body forces and con- 
sequently, the equation of state is much more satis- 
factory than with only two body force, and the sat- 
uration point becomes closer to the empirical point. 
This gives a good motivation for looking for a sim- 
ple effective interaction - or energy density functional 
- whose parameters are determined by the BHF re- 
sults. 

We have paid special attention to the effective 
mass properties in order to get constraints on the 
velocity-dependent part of the effective force. We 
find that, if a Skyrme force obeys the EBHF effective 
mass constraints then the neutron and proton effec- 
tive masses are respectively increasing and decreasing 
when the asymmetry (3 becomes larger, a behaviour 
not always obeyed by usual Skyrme forces, but con- 
sistent with empirical optical potential models[27|. 

The LNS force that we have thus obtained works 
reasonably well for describing finite nuclei in the HF 
approximation. The deviations from the data re- 
main at the level of a few percent for binding en- 
ergies and radii. If the BHF results could be im- 
proved with respect to the equilibrium point of the 
equation of state, then the HF results of finite nuclei 
would most probably become much more satisfac- 
tory. In any case, it is highly desirable to establish a 
link between microscopic many-body theories which 
are carried out in infinite systems and phenomeno- 
logical approaches for finite nuclei based on effective 
interactions. The present work is one step into this 
direction. 
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